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We study the zero-temperature crystalline structure of baby Skyrmions by applying a full-field 
numerical minimization algorithm to baby Skyrmions placed inside different parallelogramic unit- 
cells and imposing periodic boundary conditions. We find that within this setup, the minimal energy 
is obtained for the hexagonal lattice, and that in the resulting configuration the Skyrmion splits 
into quarter- Skyrmions. In particular, we find that the energy in the hexagonal case is lower than 
the one obtained on the well-studied rectangular lattice, in which splitting into half- Skyrmions is 
observed. 
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I. INTRODUCTION 

The Skyrme model [1] is a non-linear theory for pions 
in (3+1) dimensions with topological soliton solutions 
called Skyrmions. The model can be used to describe, 
with due care @, @, 3, systems of a few nucleons, and 
has also been applied to nuclear and quark matter. One 
of the most complicated aspects of the physics of hadrons 
is the behavior of the phase diagram of hadronic matter 
at finite density at low or even zero temperature. Par- 
ticularly, the properties of zero-temperature Skyrmions 
on a lattice are of interest, since the behavior of nuclear 
matter at high densities is now a focus of considerable 
interest. Within the standard zero-temperature Skyrme 
model description, a crystal of nucleons turns into a crys- 
tal of half nucleons at finite density 0, @, 0, B Q • 

To study Skyrmion crystals one imposes periodic 
boundary conditions on the Skyrme field and works 
within a unit cell (equivalently, 3-torus) T 3 ([10], p. 382). 
The first attempted construction of a crystal was by 
Klebanov [5], using a simple cubic lattice of Skyrmions 
whose symmetries maximize the attraction between near- 
est neighbors. Other symmetries were proposed which 
lead to slightly lower, but not minimal, energy crystals 
0,0. It is now understood that it is best to arrange the 
Skyrmions initially as a face-centered cubic lattice, with 
their orientations chosen symmetrically to give maximal 
attraction between all nearest neighbors |8|, [9[ . 

The Skyrme model has an analogue in (2+1) dimen- 
sions known as the baby Skyrme model, also admitting 
stable field configurations of a solitonic nature character- 
ized by integral topological charges [ll|, [l2[ • Due to its 
lower dimension, it serves as a simplification of the origi- 
nal model, but nonetheless it has a physical relevance in 
its own right, as a variant of the model arises in ferro- 
magnetic quantum Hall systems [l3|, [3, EH, Gil 03 • This 
effective theory is obtained when the excitations relative 
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to ferromagnetic quantum Hall states are described, in 
terms of a gradient expansion in the spin density, a field 
with properties analogous to the pion field in the 3D case 

The baby Skyrme model has also been studied in con- 
nection with baby Sk yrm ion lattices under various set- 
tings [H, HI E0, HI, ll and in fact, it is known that 
the baby Skyrmions also split into half- Skyrmions when 
placed inside a rectangular lattice [21]. However, to our 
knowledge, to date it hasn't been known if the rectangu- 
lar periodic boundary conditions yield the true minimal 
energy configurations over all possible lattice types. Con- 
versely, if there are other non-rectangular baby Skyrmion 
lattice configurations which have lower energy. Finding 
the answers to these questions is thus of particular im- 
portance both because of their relevance to quantum Hall 
systems in two-dimensions, and also because they may be 
used to conjecture the crystalline structure of nucleons in 
three-dimensions. 

In two dimensions there are five lattice types as given 
by the crystallographic restriction theorem [2J] , in all of 
which the fundamental unit cell is a certain type of a 
parallelogram. To find the crystalline structure of the 
baby Skyrmions, we place them inside different parallel- 
ograms with periodic boundary conditions and find the 
minimal energy configurations over all parallelograms of 
fixed area (thus keeping the Skyrmion density fixed). As 
we show later, there is a certain type of parallelogram, 
namely the hexagonal, which yields the minimal energy 
configuration. In particular, its energy is lower than the 
known 'square-cell' configurations in which the Skyrmion 
splits into half- Skyrmions. As will be pointed out later, 
the hexagonal structure revealed here is not unique to 
the present model, but also arises in other solitonic mod- 
els, such as Ginzburg-Landau vortices [25[, quantum Hall 
systems [3, EH , and even in the context of 3D Skyrmions 
[26|. The reason for this will be discussed in the conclud- 
ing section. 

In the following section we review the setup of our nu- 
merical computations, introducing a systematic way for 
the detection of the minimal energy crystalline structure 
of baby Skyrmions. In subsequent sections we outline the 
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numerical procedure through which the full field minimal 
configurations are obtained. In section HV1 we present the 
main results of our study and in section [VI a somewhat 
more analytical analysis of the problem is presented. In 
the last section we make some remarks with regards to 
future research. 



II. BABY SKYRMIONS INSIDE A 
PARALLELOGRAM 

The target manifold in the baby Skyrme model is de- 
scribed by a three-dimensional vector with the con- 
straint 0-0=1 (i.e., G S2) and its Lagrangian density 
is: 

£ = \d^-d^-^[{d^-d^) 2 (i) 

- • 3„0) • • 0»] -^ 2 (l-</> 3 ) 

consisting of a kinetic term, a 2D Skyrme term, and a 
potential term. The static solutions of the model are 
those field configurations which minimize the static en- 
ergy functional: 

E=^J^dxdy ((^0) 2 + (^0) 2 W(d x 4>xd y <f>) 2 

+2^(1-03)) (2) 

within each topological sector, where, in our setup, the 
integration is over parallelograms denoted here by A: 

A = {ai(L, 0) + a2(sLsin7, 5LCOS7) : < cei, 0L2 < 1} .(3) 

Here L is the length of one side of the parallelogram, sL 
with < s < 1 is the length of its other side and < 7 < 



7r/2 is the angle between the L sU side and the vertical to 
the 'L' side (as sketched in Fig. [I]). Each parallelogram 
is thus specified by a set {L, 5, 7} and the Skyrmion den- 
sity inside a parallelogram is p = B/(sL 2 cos 7), where B 
is the topological charge of the Skyrmion. The periodic 
boundary conditions are taken into account by identify- 
ing each of the two opposite sides of a parallelogram as 
equivalent: 

4>{x) = 4>{x + rii (L, 0) + ri2 (sL sin 7, sL cos 7)) , (4) 

with (711,712 G Z). We are interested in static finite- 
energy solutions, which in the language of differential ge- 
ometry are T2 1— > S2 maps. These are partitioned into 
homotopy sectors parameterized by an invariant integral 
topological charge £?, the degree of the map, given by: 

b = ^J a dxd V (0 • (d*0 x d y cj>)) . (5) 

The static energy E can be shown to satisfy 

E > 4tt£ (6) 
with equality possible only in the 'pure' 0(3) case (i.e., 
when both the Skyrme and the potential terms are ab- 
sent) [21.]. We note that while in the baby Skyrme model 
on R 2 with fixed boundary conditions the potential term 
is necessary to prevent the solitons from expanding in- 
definitely, in our setup it is not required, due to the pe- 
riodic boundary conditions [2l|. In this work we study 
the model both with and without the potential term. 

The problem in question can be simplified by a linear 
mapping of the parallelograms A into the unit-area two- 
torus T2 (see Fig. [1]). In the new coordinates, the energy 
functional becomes 



E = 



1 

2s cos 7 



/ dxd^( 5 2 (^0) 2 -2 5 sin 7 ^0^0 + (^0) 2 ) + ^ / dxd^/(^0x^0) 2 + — / dxdj/(l-0 3 ) -(7) 
Jt 2 ^ Jt 2 P Jt 2 



Note that the dependence of the energy on the Skyrme 
parameters k and \i and the Skyrmion density p is only 
through n 2 p and p? / 'p. 



In order to find the minimal energy configuration of 
Skyrmions over all parallelograms with fixed area (equiv- 
alently, a specified p), we scan the parallelogram parame- 
ter space {5,7} and find the parallelogram for which the 
resultant energy is minimal over the parameter space. An 
alternative approach to this problem, which is of a more 
analytical nature, has also been implemented, and is dis- 
cussed in detail in section [VI In the following section we 
describe the numerical minimization procedure. 



III. THE NUMERICAL MINIMIZATION 
PROCEDURE 



The Euler-Lagrange equations derived from the energy 
functional Q are non-linear Pi^E"s, so in general one 
must resort to numerical techniques. The minimal en- 
ergy configuration of a Skyrmion with charge B and a 
given parallelogram (defined by a set {5, 7}) with a spec- 
ified Skyrmion density p is found by a full-field relaxation 
method on a 100 x 100 rectangular grid on the two-torus 
T2, where at each point a field triplet <fi(x m ,y n ) is de- 
fined. For the calculation of the energy and charge den- 
sities, we use the procedure devised by [27], in which the 
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FIG. 1: The parameterization of the fundamental unit-cell 
parallelogram A (in black) and the two-torus T2 into which it 
is mapped (in gray). 



evaluation of these quantities is performed at the centers 
of the grid squares. Numerical derivatives are also eval- 
uated at these points; the x-derivatives are calculated by 



dx 




+ 2 n ^ 2 

/ </>(x m +i,y n ) + j)(% + l,|/n+l) 
\ 2 

0(#m, Vn) + 4>( x rn, 2/n+l) 



(8) 



normed 



normed / 

with the ^-derivatives analogously defined, and the 
"normed" subscript indicates that the averaged fields are 
normalized to one. We find this procedure to work very 
well in practice. For a more detailed discussion of the 
method, see [27j. 

The relaxation process begins by initializing the field 
triplet to a rotationally-symmetric configuration 

^initial = ( sin /( r ) cos B0 ^ sin /( r ) sin B0 ^ cos /( r )) '( 9 ) 

where the profile function f(r) is set to f(r) = 7rexp(— r) 
with r and being the usual polar coordinates. The en- 
ergy of the baby Skyrmion is then minimized by repeating 
the following steps: a point (x m ,y n ) on the grid is cho- 
sen at random, along with one of the three components 
of the field <fi(x m ,y n ). The chosen component is then 
shifted by a value 5$ chosen uniformly from the segment 
[— A^A^] where A^ = 0.1 initially. The field triplet is 
then normalized and the change in energy is calculated. 
If the energy decreases, the modification of the field is 
accepted and otherwise it is discarded. The procedure 
is repeated, while the value of A^ is gradually decreased 
throughout the procedure. This is done until no further 
decrease in energy is observed. 

To check the stability and reliability of our numerical 
prescription, we set up the minimization scheme using 
different initial configurations and an 80 x 80 grid for sev- 
eral p, s and 7 values. This was done in order to make 



sure that the final configurations are independent of the 
discretization and of the cooling scheme. As a further 
check, some of the minimizations were repeated with nu- 
merical derivatives of a higher precision, using eight field 
points for their evaluation. No apparent changes in the 
results were detected. 



IV. 



RESULTS 



Using the minimization procedure presented in the pre- 
vious section, we have found the minimal energy static 
Skyrmion configurations over all parallelograms, for var- 
ious settings: The 'pure' 0(3) case, in which both ft, the 
Skyrme parameter, and /i, the potential coupling, are set 
to zero, the Skyrme case for which only p = 0, and the 
general case for which neither the Skyrme term nor the 
potential term vanish. 

In each of these settings, the parameter space of paral- 
lelograms was scanned, while the Skyrmion density p is 
held fixed, yielding for each set of {5, 7} a minimal energy 
configuration. The choice as to how many Skyrmions are 
to be placed inside the unit cells was made after some 
preliminary testing in which Skyrmions of other charges 
(up to B = 8) were also examined. The odd-charge 
minimal-energy configurations turned out to have sub- 
stantially higher energies than even-charge ones, where 
among the latter, the charge- two Skyrmion was found 
to be the most fundamental, as it was observed that it 
serves as a 'building-block' for the higher-charge ones. In 
what follows, a summary of the results is presented. 



A. The pure 0(3) case (ft = \i — 0) 

The pure 0(3) case corresponds to setting both ft and 
p to zero. In this case, analytic solutions in terms of 
Weierstrass elliptic functions may be found [U, [22|, [23[ 
and the minimal energy configurations, in all parallelo- 
gram settings, saturate the energy bound in (j6]) giving 
E = 4ttB. Thus, comparison of our numerical results 
with the analytic solutions serves as a useful check on the 
precision of our numerical procedure. The agreement is 
to 6 significant digits. Contour plots of the charge densi- 
ties for different parallelogram settings for the charge- two 
Skyrmions are shown in Fig. O all of them of equal en- 
ergy E/8tt = 1. 



B. The Skyrme case (ft / 0, p — 0) 

As pointed out earlier, in setting p = 0, the depen- 
dence of the energy functional on the Skyrme parame- 
ter ft is only through ft 2 p, so without loss of generality 
we vary p and fix ft 2 = 0.03 throughout (this particular 
choice for ft was made for numerical convenience). Min- 
imization of the energy functional (O over all parallelo- 
grams yielded the following. For any fixed density p, the 
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functional), the Skyrmion structure is different at low 
and at high densities and a phase transition is observed. 
While at low densities the Skyrmions stand isolated from 
one another, at high densities they fuse together forming 
the quarter- Skyrmion crystal as in the Skyrme case re- 
ported above. As the density p decreases or equivalent ly 
the value of \i increases, the size of the Skyrmions be- 
comes small compared to the cell size. The exact shape 
of the lattice loses its effects and the differences in energy 
among the various lattice types become very small. This 
is illustrated in figures [4] and [5l 

Furthermore, due to the finite size of the Skyrmion, 
there is an optimal density for which the energy is min- 
imal among all densities. Figure [6] shows the contour 
plots of the charge density of the charge-two Skyrmion 
for several densities with n 2 = 0.03 and p? = 0.1. The 
energy of the Skyrmion is minimal for p « 0.14 (Fig. [4]). 



FIG. 2: Charge- two Skyrmions in the pure 0(3) case: Con- 
tour plots of the charge densities ranging from violet (low 
density) to red (high density) for various parallelogram set- 
tings, all of which saturate the energy bound E — 4ttB = 8tt. 



minimal energy was obtained for s = 1 and 7 = 7r/6. This 
set of values corresponds to the 'hexagonal' or 'equilateral 
triangular' lattice. In this configuration, any three adja- 
cent zero-energy loci (or 'holes') are the vertices of equi- 
lateral triangles, and eight distinct high-density peaks are 
observed (as shown in Fig. l3b|) . This configuration can 
thus be interpreted as the splitting of the two- Skyrmion 
into eight quarter- Skyrmions. This result turned out to 
be independent of the Skyrmion density p. In particular, 
the well-studied square-cell minimal energy configuration 
(Fig. [3aJ), in which the two-Skyrmion splits into four half- 
Skyrmions turned out to have higher energy than the 
hexagonal case. Figure [3] shows the total energies (di- 
vided by 8tt) and corresponding contour plots of charge 
densities of the hexagonal, square and other configura- 
tions (for comparison), all of them with p = 2. 

The total energy of the Skyrmions in the hexagonal 
setting turned out to be linearly proportional to the den- 
sity of the Skyrmions, reflecting the scale invar iance of 
the model (Fig. [4j). In particular, the global minimum of 
E = 4ttB = 8tt is reached when p — > 0. This is expected 
since setting p = is equivalent to setting the Skyrme 
parameter n to zero, in which case the model is effectively 
pure 0(3) and inequality (|6]) is saturated. 



C. The general case (k 7^ 0, p 7^ 0) 

The hexagonal setting turned out to be the energeti- 
cally favorable in the general case as well. Since, however, 
in this case, the Skyrmion possesses a definite size (as can 
be verified by looking at the p dependence in the energy 



V. SEMI-ANALYTICAL APPROACH 

The energy functional depends both on the Skyrme 
field 4> an d on the parallelogram parameters 7 and s. 
Formally, the minimal energy configuration over all par- 
allelograms may be obtained by functional differentia- 
tion with respect to 4> an d regular differentiation with 
respect to 7 and s. However, since the resulting equa- 
tions are very complicated, a direct numerical solution is 
quite hard. Nonetheless, some analytical progress may 
be achieved in the following way. As a first step, we dif- 
ferentiate the energy functional ([7j) only with respect to 




(a) a = 1, 7 = 0, 
and E/8tt = 1.446 



Mill 

(c) s = 0.51, 7 = 0, and 
E/8tt = 1.587 




(b) s = 1,7 = 7r/6, and 
E/8tt = 1.433 



(d) s = 1, 7 = 7r/4, and 
E/8tt = 1.454 



FIG. 3: Charge- two Skyrmions in the Skyrme case with 
k 2 = 0.03 and p — 2: Contour plots of the charge densities 
for the hexagonal, square and other fundamental cells rang- 
ing from violet (low density) to red (high density). As the 
captions of the individual subfigures indicate, the hexagonal 
configuration has the lowest energy. 
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FIG. 4: Total energy E (divided by 8tt) of the charge-two 
Skyrmion in the hexagonal lattice (□ - Skyrme case and 
- general case) and in the square lattice ( ■ - Skyrme case 
and ♦ - general case) as function of the Skyrmion density (in 



the Skyrme case, k — 0.03 and in the general case k 



0.03 



and /j 2 = 0.1). Note the existence of an optimal density (at 
p « 0.14) in the general case, for which the energy attains a 
global minimum. Figure (b) is an enlargement of the lower 
left corner of figure (a). 



7 and s: 
dE 



1 



c^7 2s cos 2 7 
dE _ 1 
ds 2s 2 cos 7 



(sm-f(£ yy + s 2 £ xx ) - 2s£ xy ) = , 



(£ yy -s 2 £ xx ) = 0, 



(10) 



where £ l i = f T dxdy(di<f)'dj<f)) andi,j G {x,y}. Solving 
these equations for 7 and s yields 



sin 7 



£yy 

£xx ' 
gxy 



(11) 



\/£ xx £yy ' 



Substituting these expressions into the energy functional 
(0), we arrive at a 'reduced' functional 



E 



2B 



pot • 



(12) 



where £ s ^ = J T2 dxdy (d x 4> x dy(f)) 2 is the Skyrme energy 
and £pot = Jj 2 dxdy (1 — 03) is the potential energy. 
Now that both 7 and s are eliminated from the resultant 
expression, and the conditions for their optimization are 
built into the functional, the numerical minimization is 
carried out. We note here, however, that the procedure 
presented above should be treated with caution. This is 
since the 'minimization' conditions for 7 and s in ([TTj) are 
in fact only extremum conditions, and may turn out to be 
maximum or saddle-point conditions. Hence, it is impor- 
tant to confirm any results obtained using this method 
by comparing them with corresponding results obtained 
from the method described in the previous section. 

It is therefore reassuring that numerical minimization 
of the reduced functional ([T2]) gives sin 7 = 0.498 (7 « 
7r/6) and s = 1 (both for the Skyrme case and the general 
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FIG. 5: Energy difference (in %) between the hexagonal 
setting and the square setting, in the Skyrme case (□) and in 
the general case (♦), as function of the Skyrmion density. As 
the density decreases, the energy difference between the two 
lattice types is reduced. 
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(a) p = 1, E/8n = 1.229 (b) p = 0.5, E/8ir = 1.128 




(c) p = 0.25, £/8tt - 1.087 (d) p = 0.14, E/8tt = 1.08 




(e) p = 0.0625, £/8tt = 1.082 (f) p = 0.03125, E/8ir = 1.084 



FIG. 6: Charge- two Skyrmions in the general case with 
/j 2 = 0.1 and k 2 = 0.03: Contour plots of the charge den- 
sities of the minimal-energy configurations in the hexagonal 
setting for different densities. Here, the energetically most fa- 
vorable density is p — 0.14. The plot colors range from violet 
(low density) to red (high density). 

case), confirming the results presented in the previous 
section. 

In the general (/i =^ 0) case, the energy functional 
([T2]) may be further differentiated with respect to the 
Skyrmion density p to obtain the optimal density for 
which the Skyrmion energy is minimal. Differentiating 
with respect to p, and substituting the obtained expres- 
sion into the energy functional, results in the functional 

Numerical minimization of the above expression for n 2 = 
0.03 and various p values (0.1 < p 2 < 10) yielded the 
hexagonal setting as in the Skyrme case. In particular, 
for p 2 — 0.1 the optimal density turned out to be p ~ 
0.14, in accord with results presented in sub-section llV CI 

VI. SUMMARY AND FURTHER REMARKS 

In this paper, we have studied the crystalline struc- 
ture of baby Skyrmions in two dimensions by finding the 
minimal energy configurations of baby Skyrmions placed 
inside parallelogramic fundamental unit cells and impos- 
ing periodic boundary conditions. In the pure 0(3) case 
(where both the Skyrme and potential terms are absent), 



we verified that there are no favorable lattices in which 
the Skyrmions order themselves, as all parallelogram set- 
tings yielded the same minimal energy, saturating the 
minimal energy bound given by inequality ([6|). 

In the Skyrme case, without the potential term, the 
results are different. For any fixed Skyrmion density, the 
parallelogram for which the Skyrmion energy is minimal 
turns out to be the hexagonal lattice (for which s = 1 and 
7 = 7r/6). In particular, the hexagonal lattice has lower 
energy than the four half- Skyrmions configuration on the 
rectangular lattice. For example, at p = 2, E — 1.433 for 
the hexagonal lattice vs. E = 1.446 for the rectangular 
lattice. 

The hexagonal setting turns out to be the energetically 
most favorable also in the general case with both the 
Skyrme term and the potential term. In this case, how- 
ever, the model is not scale invariant and the Skyrmion 
has a definite size. This results in the existence of a 
phase transition as a function of density and the appear- 
ance of an optimal density, for which the total energy 
of the Skyrmion is minimal. This is analogous to the 
Skyrmion behavior in the 3D Skyrme model, in which 
the Skyrmions also possess a definite size [5[. 

As pointed out in the Introduction, the special role of 
the hexagonal lattice revealed here is not unique to the 
baby Skyrme model, but in fact arises in other solitonic 
models. In the context of Skyrme models, the existence of 
a hexagonal two dimensional structure of 3D Skyrmions 
has been found by [26[ , where it has already been noted 
that energetically, the optimum infinite planar structure 
of 3D Skyrmions is the hexagonal lattice, which resem- 
bles the structure of a graphite sheet, the most stable 
form of carbon thermodynamically ([lQ[, p. 384). Other 
examples in which the hexagonal structure is revealed are 
Ginzburg-Landau vortices which are known to have lower 
energy in a hexagonal configuration than in a square lat- 
tice configuration [25[. Thus, it should not come as a 
surprise that the hexagonal structure is found to be the 
most favorable in the baby Skyrme model. 

As briefly noted in the Introduction, a certain type 
of baby Skyrmions also arise in quantum Hall systems as 
low-energy excitations of the ground state near ferromag- 
netic filling factors (notably 1 and 1/3) [141 ]. It has been 
pointed out that this state contains a finite density of 
Skyrmions [28[, and in fact the hexagonal configuration 
has been suggested as a candidate for their lattice struc- 
ture [15[ . Our results may therefore serve as a supporting 
evidence in that direction, although a more detailed anal- 
ysis is in order. 

Our results also raise some interesting questions. The 
first is how the dynamical properties of baby Skyrmions 
on the hexagonal lattice differ from their behavior in the 
usual rectangular lattice. Another question has to do 
with their behavior in non-zero-temperature settings. 

One may also wonder whether and how these results 
can be generalized to the 3D case, once a systematic 
study like the one reported here is conducted. Is the 
face-centered cubic lattice indeed the minimal energy 
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crystalline structure of 3D Skyrmions among all paral- 
lelepiped lattices? If not, what would the minimal en- 
ergy structure be? and how would these results change 
when a mass term is present? We hope to answer these 
questions in future research. 
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